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Abstract 



"Constants of Nature" and cosmological parameters may in fact be vari- 
ables related to some slowly- varying fields. In models of eternal inflation, such 
fields will take different values in different parts of the universe. Here I show 
how one can assign probabilities to values of the "constants" measured by 
a typical observer. This method does not suffer from ambiguities previously 
discussed in the literature. 
One of the striking aspects of the inflationary cosmology is that, generically, inflation 
never ends p|. The evolution of the infiaton field is influenced by quantum fluctuation, and 
as a result thermalization does not occur simultaneously in different parts of the universe. 
The dynamics of can be pictured as a random walk superimposed on the classical slow roll 
and can be described by a Fokker-Planck equation IQ-Q. On small scales, the fluctuations 
in the thermalization time give rise to a spectrum of small density fluctuations, but on large 
scales they make the universe extremely inhomogeneous. In most of the models one flnds 
that at any time there are parts of the universe that are still inflating and that the total 
volume of the inflating regions is growing with time. This picture is often referred to as 
"stochastic inflation" or "eternal inflation". 

Eternally inflating universes may contain thermalized regions characterized by different 
values of the constants of Nature and of the cosmological parameters. For example, the in- 
fiaton potential V{(f)) may have several minima corresponding to different low-energy physics 
or to different values of the cosmological constant. A more interesting possibility is that the 
"constants" are related to some slowly-varying fields and take values in a continuous range. 
Examples are the effective gravitational constant in a Brans-Dicke-type theory and the 
density parameter f2 in some models of open inflation 0. Just like the inflaton fleld 0, the 
flelds XjU = 1, •••,'"') associated with the "constants" are subject to quantum fluctuations 
during inflation, and different regions of the universe thermalize with different values of Xj- 
An intriguing question is whether or not we can predict which values of the "constants" we 
are most likely to observe. 

An eternally inflating universe is inhabited by numerous civilizations that will measure 
different values of Xj- (Fo^^ simplicity, I do not distinguish between the flelds Xj ^ind the 
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associated constants). We can define the probability V{x)d"'X ^^i Xj to be in the intervals dxj 
as being proportional to the number of civilizations which will measure Xj in that interval [0 . 
Assuming that we are a typical civilization, we can expect to observe Xj near the maximum 
oiVix) i- 

An immediate objection to this approach is that we are ignorant about the origin of 
life, let alone intelligence, and therefore the number of civilizations cannot be calculated. 
However, the approach can still be used to find the probability distribution for parameters 
which do not affect the physical processes involved in the evolution of life. The cosmolog- 
ical constant A, the density parameter Q and the amplitude of density fluctuations Q are 
examples of such parameters. We shall assume that our fields Xj belong to this category. 
The probability for a civilization to evolve on a suitable planet is then independent of Xj, 
and instead of the number of civilizations we can use the number of habitable planets or, as 
a rough approximation, the number of galaxies. 

If this philosophy is adopted, the problem of calculating V{x) can be split into two 
parts. The probability for us to observe certain values of Xj is proportional to the volume 
of the regions where Xj take specified values and to the density of galaxies in those regions. 
It is convenient to consider comoving regions and measure their volume at the time of 
thermalization. Then we can write 



^(X)0cz/(X)^*(X)- (1) 

Here, is proportional to the volume of thermalized regions where Xj take values in 

the intervals dxj, and is the number of galaxies that form per unit thermalized volume. 
The calculation of is a standard astrophysical problem, which is completely unrelated 
to the calculation of the volume factor Our focus in this paper will be on the volume 

factor. 

The problem one encounters in using Eq.(l) is that in an eternally inflating universe the 
thermalized volume is infinite, even for a region of a finite comoving size. One can deal with 
this problem by simply introducing a time cutoff and including only regions that thermalized 
prior to some moment of time tc- One finds, however, that the resulting probability distri- 
bution is extremely sensitive to the choice of the time variable t 0]. This gauge- dependence 
casts doubt on any conclusion reached using this approach. 

An alternative procedure, suggested in is to introduce a x-dependent cutoff at 

the time tc{x), when all but a small fraction e of the comoving volume destined to thermalize 
with Xj in the intervals dxj has thermalized. The limit e — > is taken after calculating the 
probabilities. It was shown in |TD|,|TT| that the resulting probability distribution is essentially 



insensitive to the choice of time parametrization. However, the same problem appears 
in a different guise. Linde and Mezhlumian [|T^] have found a family of gauge-invariant 
cutoff procedures which includes the e-procedure described above. This indicates that the 
invariance requirement alone is not sufficient to define the probabilities uniquely. 

Another source of uncertainty is associated with large values of the infiaton potential 
V{(f)). Usually, (p is assumed to vary in an infinite range, with V{(f)) — *• oo at — >• ±cxo. 
Higher values of V{(f)) correspond to higher rates of expansion, and thus the regions with 
Planckian and super-Planckian values of V{(f)) expand at the highest rate. The effect of 
these regions on the dynamics of stochastic inflation is therefore very important |T3I]. Since 



the physics at Planckian energy density is completely unknown, this results in a significant 
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uncertainty 



In the face of these difficulties, serious doubts have been expressed that a meaningful 
definition of probabilities in an eternally infiating universe is even in principle possible 
Pj9|Jl^. However, these pessimistic conclusions may have been premature. In this paper, I 
am going to argue that unambiguous probabilities can in fact be defined in a wide class of 
models of eternal infiation. 

The factor in Eq. (|I|) is the probability distribution of the fields Xj the ther- 

malization hypersurface which separates infiating and thermalized spacetime regions. It 
is a 3-dimensional spacelike surface which plays the role of the big bang for the thermal- 
ized regions. In the case of several discrete vacua, consists of a number of disconnected 
pieces, each connected component corresponding to thermalization into a single vacuum 
1^ . Each connected piece of has an infinite volume. (The situation here is somewhat 



similar to open infiation, where thermalized regions inside the bubbles have the form of 
infinite, open Robertson- Walker universes). In order to determine the relative probability of 
different vacua, one has to compare the infinite volumes of the corresponding components 
of S*, which is an inherently ambiguous task. This ambiguity is the source of the problems 
encountered in Refs. i,|g-|T2|. 

Now, the key observation is that the situation may be greatly improved in the case of 
continuous fields Xj varying in a finite range, < Xj ^j- In this case, each connected 
part of S^, is still infinite, but now different parts are not characterized by different values 
of Xj- On the contrary, Xj run through all the range of their values on each connected part. 
Since the inflationary dynamics of the flelds Xj has a stochastic nature, the distributions of 
Xj on different components of E^, should be statistically equivalent. It is therefore sufficient 
to consider a single connected component. Moreover, since Xj have a flnite range, they will 
run through all of their values many times on any sufficiently large part of S*. Hence, there 
is no need to deal with infinite hypersurfaces. The probability distribution can be 

determined by examining a large finite piece of S*. 

The above argument is an adaptation, to the case of stochastic infiation, of the argument 
given earlier by Garriga, Tanaka and Vilenkin [15| for models of open infiation with a variable 
Q. In such models, each bubble contains regions with all possible values of fl, and we argued 
that the probability distribution for Q can be determined by considering a single bubble. 

The argument can also be extended to fields with an infinite range of variation, pro- 
vided that the probability distributions of Xj are concentrated within a finite range, with a 
negligible probability of finding Xj very far away from that range. 

Models with a finite range of Xj are not difficult to construct. For example, Xj could 
play the role of angular variables, with the infiaton field being the radial variable in the 
field space. Another attractive possibility is a a-model-type theory in which both Xj and 
(j) take values on a compact manifold. In such models the potential is necessarily bounded 
and problems with super-Planckian densities are easily avoided. 

The gauge- dependence of the probability distributions obtained using a constant-time 
cutoff can be understood as follows. Suppose that the time coordinate t has been chosen to 
be the proper time, t = r, and that regions with < xi < Ai/2 take a shorter time r to 
thermalize than regions with Ai/2 < xi < ^i- Then, in the vicinity of the cutoff (t = tc), 
the surface T,^:(tc) will tend to include regions with xi < ■^i/S and exclude regions with 
Xi > Ai/2. Suppose now that, despite the shorter roll-down time, regions with xi < ^i/2 
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expand by a larger factor than regions with xi > ^i/2. (This can be arranged with a 
suitable choice of the potential V{(f), x))- Then, if we choose the scale factor to be the time 
coordinate, t = a, the surface S^(tc) will tend to include regions with xi > ^i/2 and exclude 
those with xi < ^i/2. Thus, the cutoff at a fixed time tc is biased with respect to different 
types of regions, depending on the choice of the time coordinate. One might think that the 
probability distribution of x on a hypersurface should be insensitive to variation of the 
boundary of S^,, provided that the volume of is sufficiently large. However, in the case 
of a fixed-time cutoff, S^, is a multi-component hypersurface whose volume is dominated 
by small, newly formed thermalized regions. It is therefore no surprise that the resulting 
distribution is sensitive to the choice of the time variable. 

Thus, what we need is an unbiasly selected portion of S^,. Such an unbiased cutoff 
is hard to implement in the Fokker- Planck formalism which deals with probability 
distributions on equal-time surfaces. A more promising approach to the calculation of 
is to use a numerical simulation of stochastic inflation. Simulations of this sort have been 



introduced in Ref. |jT6[ and have been further developed in Ref. 0]. In the present paper, I 
will not attempt to pursue this numerical approach. Instead, I am going to argue that in a 
wide class of models the probabilities can be calculated analytically. 

Suppose the maximum of the potential V{(j), x) is at = 0. We shall first consider models 
in which V{(j), x) is very symmetric near 

the top, so that it is essentially independent of x i^i the range < < 0i, but may 
have significant x-dependence elsewhere. We shall assume that 0i is in the slow-roll range, 
0g -C 01 < 0*(x)) where 0^ is the boundary between the quantum diffusion and deterministic 
slow roll regimes, and 0*(x) is the value corresponding to the end of inflation. Let us consider 
a constant- hypersurface Eq: = 0o with 0g -C 0o < 0i. Since 0o is in the slow-roll range, 
this hypersurface is (almost everywhere) spacelike |T^. Moreover, all values of in the past 



of So belong to the range < < 0o where the potential is independent of x- Therefore, we 
expect the distribution of Xj on Sq to be fiat, 'Po(Xi) = const, with all values of Xj equally 
probable |T^. 

The hypersurface S^. is a spacelike hypersurface in the future of Sq. The probability 
distribution for Xj on is related to that on Sq by 

n(x) oc Po(x^°^) exp[3iV(x(°))]det(ax(°V5x). (2) 

Here, xf\x) the values on Eq that evolve, along the slow-roll trajectory, to the values Xj 
on S*, N(x^^^) is the number of inflationary e-folding along that trajectory, exp(3A'") is the 
volume expansion factor, and det((9x''°V'9x) is the Jacobian transforming from variables xf^ 
to Xj- Now, Vo{x^'^^) = const, and the functions X^^Hx) ^^nd N(x^^^) can be straightforwardly 
determined from a given potential V{(j),x)- Hence, Eq. (0) solves the problem of finding 
the probability distribution for Xj on S^,. This equation is significantly simplified in models 
where Xj do not change much during the slow roll, so that xf^ ~ Xj- such models, 

P,(x)ocexp[3iV(x)]. (3) 

Eq.(2) disregards quantum fluctuations of Xj and in the slow-roll regime at > 0o. 
For Xj this is justified, provided that the average cumulative fluctuation is small compared 
to Aj. This condition is not difficult to satisfy. Small quantum fluctuations of are not 
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correlated with Xj and their effect vanishes to hnear order in 5(j). However, quadratic and 
higher-order terms will give some corrections to Eq.(2). One may also be worried about 
the effect of very large fluctuations which bring back to the quantum diffusion regime 
Although extremely rare, such fluctuations are likely to be present for a sufficiently 
large piece of Sq. If the probability of large fluctuations were independent of Xj: then small 
patches of Sq which have such fluctuations in their future could be left out, without any 
effect on the distribution (2). If there is some x-dependence, its effect will be to reduce 
for those values of Xj which give a higher probability of large backwards jumps. I expect 
this effect to be very small, but more work is needed to obtain a quantitative estimate. 



Large fluctuations back to the quantum regime may occur even after thermalization |]18 
Then each thermalization surface may have an inflnite number of other thermalization 
surfaces in its future. However, Eq.(2) can still be used to derive the probability distribution 
for Xj on any of these surfaces. 

As an illustrative example, let us consider a two-fleld model with a potential 

V{(f),x) = Vocos^{a(f))[l + Xsm\a(j))cos{(3x)]- (4) 

Here, |A| < 1, a = vr/A^ and the range of is < < A^, with the values = and 
= A^ identifled. Similarly, (3 = vr/A^ with x = and x = identifled. We shall assume 
that j3 ^ a ^ 1, that is, A^ > A,^ ^ 1. (I use Planck units throughout the paper). The top 
of the potential is at = 0, and the true vacuum is at = The fleld x is massless in 

the true vacuum, while the mass of is x-dependent, 

ml{x) = 2Voa^[l + Xcos{f3x)]. (5) 

The boundary 0^ between the quantum diffusion and slow roll regimes is determined by 
the condition IdH/dcp] ~ H^, where -ff(0, x) = [§^^(0) x)/3]^''^ is the inflationary expansion 
rate. We shall assume that Vq <^ a^; then 0^ ~ Vq^"^ /a'^ ^ A^. Since a ^ 1, the problem 
of super-Planckian energies does not arise in this case. 

The potential (^) is essentially independent of x for -C A^, but the number of e-foldings 
and the amplitude of density fluctuations Sp/p are both x-dependent. For a <C 1, A <^ 1 
we flnd 

and i 

5p/p{x) ~ 200m^(x). (7) 

Here, A"o is the number of e-foldings (from a specifled reference point 0o) in the X" 
independent case (A = 0), and I have used sin(a0*) ^ 1, which is valid for a < 1 (since 
0^, ^ ^ ~ |)- The classical and quantum variations of x during the relevant part of the slow 

roll are, respectively, dxd ~ ('^/47i")(/5/a)^ ^ 1 and 8xq/ ^x ^ ^o^'^f^l^ 1; and can 
be neglected. This justifles the x = const approximation in Eq. (H) and the use of Eq. (^ 
for the probabihty distribution, 

'P*(x) oc exp[67rAa"^ cos(/5x)]- (8) 
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The corresponding distribution for 5p/ p can easily be written using Eqs. (||), (0). 

As a second example, we consider the spectrum of density perturbations in the standard 
model of inflation with a single field 0. The perturbations are determined by quantum 
fluctuations 50 which can be regarded as random Gaussian variables with a dispersion 
a = H/27r. A perturbation is produced on each comoving scale at the time when that scale 
crosses the horizon and has a gauge-invariant amplitude [l^ 



6p/p ~ 8nH5(l)/H', (9) 
where H' = dH/d(j). Averaging over the Gaussian distribution gives the standard result 

{5p/p)rms--^Hy\H'\. 

Fluctuations of on different length scales are statistically independent and can be 
treated separately. We can therefore concentrate on a single scale corresponding to some 
value = 00, disregarding all the rest. The fluctuation 6(f) will take different values on 
different parts of the hypersurface Eq : = 0o, so 50 plays the role of the field x of our 
previous example. Linde et. al. have studied the probability distribution for the observed 
values of 50 using a constant proper time cutoff and arrived at a surprising conclusion that 
typical fluctuations can be much greater than the rms value H/2t^. We shall see, however, 
that this large effect disappears in the present approach. 

The probability distribution for 50 on can be found directly from Eq. @ with 'Po((^0) 
being the initial Gaussian distribution on Eg. When the field undergoes a quantum jump 
50, the number of e-foldings necessary to complete the classical rollover to 0^, is changed by 
5N = AnScpHo/ H'q. Hence, we obtain the distribution 



oc exp 



-^f50-MV 



(10) 



which describes fluctuations with a nonzero mean value, < 50 >= 3Hq/ttH'q. The inflaton 
tends to fluctuate in the direction opposite to the classical roll, because backward fluctuations 
prolong inflation and increase the volume factor. This effect, however, is hopelessly small: 
< 50 > /cr = 6Hq/H'q ~ {Sp/p)rms 1- We thus see that anthropic considerations do not 



substantially modify the standard calculations of the density perturbation spectrum pT |. 

In conclusion, the approach described in this paper allows one to assign unambiguous 
probabilities to continuously-varying "constants" measured by a typical observer. In mod- 
els where the inflaton potential is nearly independent of the "constants" in the quantum 
diffusion range, the probabilities can be calculated analytically. Otherwise, a numerical 
simulation is needed to determine the distribution Po(x*^°^) iii Eq. (||). 

I am grateful to Jaume Garriga for useful comments on the manuscript. This work was 
supported in part by the National Science Foundation. 



6 



REFERENCES 



[1] For a review of inflation, see, e.g., A.D. Linde, Particle Physics and Inflationary Cos- 
mology (Harwood Academic, Chur, Switzerland, 1990); K.A. Olive, Phys. Rep. 190, 
307 (1990). 

[2] A. Vilenkin, Phys. Rev. D27, 2848 (1983). 

[3] A. A. Starobinsky, in Lecture Notes in Physics Vol. 246 (Springer, Heidelberg, 1986). 
[4] A. D. Linde, D. A. Linde, and A. Mezhlumian, Phys. Rev. D49, 1783 (1994). 
[5] J. Garcia-BeUido, A. D. Linde, and D. A. Linde, Phys. Rev. D50, 730 (1994). 
[6] A. D. Linde and A. Mezhlumian, Phys. Rev. D52, 6789 (1995). 
[7] A. Vilenkin, Phys. Rev. Lett. 74, 846 (1995). 

[8] This version of the anthropic principle, which I called the principle of mediocrity [0, 
is an extension of the Copernican principle. The Copernican principle asserts that our 
position in space is not special, while the principle of mediocrity asserts that the values 
of the cosmological parameters we are going to measure are not special either. A similar 
approach was used by B. Carter [unpublished], J. Leslie [Mind 101.403, 521 (1992)] and 
J.R. Gott [Nature 363, 315 (1993)] to estimate the expected lifetime of our species. 
Related ideas have also been discussed by A. Albrecht [in Lecture Notes in Physics, 
Springer- Verlag, Berlin, 1995] and by A.D. Linde et. al. 
[9] J. Garcia-Belhdo, A.D. Linde and D.A. Linde, Phys. Rev. D50, 730 (1994); J. Garcia- 
BeUido and A.D. Linde, Phys. Rev. D51, 429 (1995) 

[10] A. Vilenkin, Phys. Rev. D52, 3365 (1995). 

[11] S. Winitzki and A. Vilenkin, Phys. Rev. D53, 4298 (1996). 

[12] A. D. Linde and A. Mezhlumian, Phys. Rev. D53, 4267 (1996). 

[13] A. D. Linde, Phys. Lett. B175, 395 (1986). 

[14] I assume that the discrete vacua cannot be connected by non-infiating domain walls. 
See A.D. Linde, Phys. Lett. B 327, 208 (1994); A. Vilenkin, Phys. Rev. Lett. 72, 3137 
(1994). 

[15] J. Garriga, T. Tanaka and A. Vilenkin, |astro-ph / 9803268] . 
[16] M. Aryal and A. Vilenkin, Phys. Lett. B199, 351 (1987). 

[17] If the fields Xj ^ire described by a a-model with a non-trivial metric in the x-space, so 
that the kinetic terms in the Lagrangian are of the form Kij{x)d^Xid^Xjj then Vo{x) oc 
[det{Kij)Y/\ In this paper I assume for simplicity that Kij = {l/2)6ij. 

[18] J. Garriga and A. Vilenkin, Phys. Rev. D57, 2230 (1998). 

[19] For a review of density fluctuations in inflationary scenarios, see, e. g., V. F. Mukhanov, 

H. A. Feldman, and R. H. Brandenberger, Phys. Rep. 215, 203 (1992). 
[20] A. D. Linde, D. A. Linde, and A. Mezhlumian, Phys. Rett. B345, 203 (1995); Phys. 

Rev. D54, 2504 (1996). 
[21] It should be noted that our results (^ and ([lOD are in agreement with those obtained 

using the e-prescription of Ref. |Ty]. However, this agreement may not persist in the 

general case. 



7 



